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Abstract 

Hydrodynamic derivation of the entrainment of matter induced by a surface elastic wave 
propagating along the flexible vacuum-matter interface is conducted by considering the nonlinear 
coupling between the interface and the rarefaction effect. The critical reflux values associated 
with the product of the second-order (unit) body forcing and the Reynolds number (representing 
the viscous dissipations) decrease as the Knudsen number (representing the rarefaction measure) 
increases from zero to 0.1. We obtained the critical bounds for matter-freezed or zero-volume- 
flow-rate states corresponding to specific Reynolds numbers (ratio of wave inertia and viscous 
dissipation effects) and wave numbers which might be linked to the dissipative evolution of the 
Universe. Our results also show that for certain time-averaged evolution of the matter (gas) there 
might be existence of negative pressure. 
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I. INTRODUCTION 



The mean cosmic density of dark matter (plus baryons) is now pinned down to be only 
ca. 30% of the so-called critical density corresponding to a 'flat'-Universe. However, other 
recent evidence — microwave background anisotropies, complemented by data on distant su- 
pernovae — reveals that our Universe actually is 'flat', but that its dominant ingredient (ca. 
70% of the total mass energy) is something quite unexpected: 'dark energy' pervading all 
space, with negative pressure. We do know that this material is very dark and that it 
dominates the internal kinematics, clustering properties and motions of galactic systems. 
Dark matter is commonly associated to weakly interacting particles (WIMPs), and can be 
described as a fluid with vanishing pressure. It plays a crucial role in the formation and 
evolution of structure in the universe and it is unlikely that galaxies could have formed 
without its presence [1-2]. Analysis of cosmological mixed dark matter models in spatially 
flat Friedmann Universe with zero A term have been presented before. For example, we 
can start from the Einstein action describing the gravitational forces in the presence of the 
cosmological constant [3] 



2k N 

where S m is the contribution of the matter and radiation, k 2 N = 8tiGn = 8n/Mp, is the 
Newton's constant and Mp = 1.22 x 10 19 GeV/c 2 is the Planck mass. The set of equations 
governing the evolution of the universe is completed by the Friedman equations for the scale 
factor (a(t)) 

(«)> = !!«!*(, + A)* 



a )2 = _4^ _ 
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where k = is for the flat- Universe, p and p is the density and pressure [4-5]. A large ma- 
jority of dark energy models describes dark energy in terms of the equation of state (EOS) 
Pd = uj pa, where uj is the parameter of the EOS, while pd and pd denote the pressure and the 
energy density of dark energy, respectively. The value uj = — 1 is characteristic of the cos- 
mological constant, while the dynamical models of dark energy generally have uj > —1. The 
case of the growing cosmological term A and its implications for the asymptotic expansion of 
the universe and the destiny of the bound systems have been studied in [5] using above sys- 
tem of equations. Their results showed that even for very slow growth of A (which satisfies 



all the conditions on the variation of Gn), in the distant future the gravitationally bound 
systems become unbound, while the nongravitationally bound systems remain bound. 
Meanwhile, it is convenient to express the mean densities pi of various quantities in the 
Universe in terms of their fractions relative to the critical density: fij = Pi/p cr it- The theory 
of cosmological inflation strongly suggests that the total density should be very close to the 
critical one (Qtat ~ 1), an d this is supported by the available data on the cosmic microwave 
background (CMB) radiation. The fluctuations observed in the CMB at a level ca. 1CT 5 
in amplitude exhibit a peak at a partial wave / ~ 200, as would be produced by acoustic 
oscillations in a flat Universe with Q tat ~ 1. At lower partial waves, / 3> 200, the CMB 
fluctuations are believed to be dominated by the Sachs- Wolfe effect due to the gravitational 
potential, and more acoustic oscillations are expected at I > 200, whose relative heights 
depend on the baryon density fV At even larger values, I > 1000, these oscillations should 
be progressively damped away [1-2,3,6]. 

Influential only over the largest of scales-the cosmological horizon-is the outermost species of 
invisible matter: the vacuum energy (also known by such names as dark energy, quintessence, 
x-matter, the zero-point field, and the cosmological constant A). If there is no exchange of 
energy between vacuum and matter components, the requirement of general covariance im- 
plies the time dependence of the gravitational constant G. Thus, it is interesting to look at 
the interacting behavior between the vacuum (energy) and the matter from the macroscopic 
point of view. One related issue, say, is about the dissipative matter of the flat Universe 
immersed in vacua [7] and the other one is the macroscopic Casimir effect with the deformed 
boundaries [8]. 

Theoretical (using the Boltzmann equation) and experimental studies of interphase nonlocal 
transport phenomena which appear as a result of a different type of nonequilibrium repre- 
senting propagation of a surface elastic wave have been performed since late 1980s [9-10]. 
These are relevant to rarefied gases (RG) flowing along deformable elastic slabs with the 
dominated parameter being the Knudsen number (Kn = mean-free-path/L^, mean-free-path 
(mfp) is the mean free path of the gas, Ld is proportional to the distance between two slabs) 
[11-13]. The role of the Knudsen number is similar to that of the Navier slip parameter 
N s [14]; here, N s = pS/d is the dimensionless Navier slip parameter; S is a proportionality 
constant as u s = St, t : the shear stress of the bulk velocity; u s : the dimensional slip 
velocity; for a no-slip case, 5 = 0, but for a no-stress condition. S — oo, p is the fluid 



3 



viscosity, d is one half of the distance between upper and lower slabs). 
Here, the transport driven by the wavy elastic vacuum-matter interface will be presented. 
The flat-Universe is presumed and the corresponding matter is immersed in vacua with the 
interface being flat-plane like. We adopt the macroscopic or hydrodynamical approach and 
simplify the original system of equations (related to the momentum and mass transport) 
to one single higher-order quasi-linear partial differential equation in terms of the unknown 
stream function. In this study, we shall assume that the Mach number Ma <C 1, and 
the governing equations are the incompressible Navier-Stokes equations which are associ- 
ated with the relaxed slip velocity boundary conditions along the interfaces [11-14]. We 
then introduce the perturbation technique so that we can solve the related boundary value 
problem approximately. To consider the originally quiescent gas for simplicity, due to the 
difficulty in solving a fourth-order quasi-linear complex ordinary differential equation (when 
the wavy boundary condition are imposed), we can finally get an analytically perturbed 
solution and calculate those physical quantities we have interests, like, time-averaged trans- 
port or entrainment, perturbed velocity functions, critical unit body forcing corresponding 
to the freezed or zero-volume-flow-rate states. These results might be closely linked to the 
vacuum-matter interactions (say, macroscopic Casimir effects) and the evolution of the Uni- 
verse (as mentioned above : the critical density [1-2]). Our results also show that for certain 
time-averaged evolution of the matter (the maximum speed of the matter (gas) appears at 
the center-line) there might be existence of negative-pressure states. 



II. FORMULATIONS 



We consider a two-dimensional matter-region of uniform thickness which is approximated 
by a homogeneous rarefied gas (Newtonian viscous fluid). The equation of motion is 

(A L + /i)graddivu + /iVu + pp = p— , (1) 

where and /x are Lame constants, u is the displacement field (vector), p is the mass density 
and p is the body force for unit mass. The Navier-Stokes equations, valid for Newtonian 
fluids (both gases and liquids), has been a mixture of continuum fluid mechanics ever since 
1845 following the seemingly definitive work of Stokes and others [15], who proposed the 
following rheological constitutive expression for the fluid deviatoric or viscous stress (tensor) 
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T : T = 2/iVu + A L IV ■ u. 

The flat-plane boundaries of this matter-region or the vacuum-matter interfaces are rather 
flexible and presumed to be elastic, on which are imposed traveling sinusoidal waves of 
small amplitude a (possibly due to vacuum fluctuations). The vertical displacements of the 
upper and lower interfaces (y — h and —h) are thus presumed to be i] and — rj, respectively, 
where 77 = acos[27r(a; — ct)/X], A is the wave length, and c the wave speed, x and y are 
Cartesian coordinates, with x measured in the direction of wave propagation and y measured 
in the direction normal to the mean position of the vacuum-matter interfaces. It would be 
expedient to simplify these equations by introducing dimensionless variables. We have a 
characteristic velocity c and three characteristic lengths a, A, and h. The following variables 
based on c and h could thus be introduced : 

, x , y , u , v , 77 , -0 ./ ct , p 

X = T, y = T, U = -, V = -, 77 = -, 1p = — , t = — , p = — , 

h h c c h ch h pc 1 

where ip is the dimensional stream function, u and v are the velocities along the x- and 
y-directions; p is the density, p (its gradient) is related to the (unit) body forcing. The 
primes could be dropped in the following. The amplitude ratio e, the wave number a, and 
the Reynolds number (ratio of wave inertia and viscous dissipation effects) Re are defined 
by 

a 2nh ch 

e — — , a — — — , Re = — . 
h A v 

We shall seek a solution in the form of a series in the parameter e : 

i> = + e^i + e 2 V> 2 H , 

dp dp dp 2 dp 

with u = dip/dy, v = —dip/dx. The 2D (x- and y-) momentum equations and the equation 
of continuity could be in terms of the stream function ip if the p-term (the specific body 
force density, assumed to be conservative and hence expressed as the gradient of a time- 
independent potential energy function) is eliminated. The final governing equation is 

and subscripts indicate the partial differentiation. Thus, we have 

d 1 

— V 2 ^o + ^o y V 2 Vox - V'oxV 2 ^ = ^VVo, (3) 



d 1 
d 

and other higher order terms. The (matter) gas is subjected to boundary conditions imposed 
by the symmetric motion of the vacuum-matter interfaces and the non-zero slip velocity : 
u = =1= Kn du/dy [11-13], v = ±dr]/dt at y = ±(1 + 1]), here Kn=mfp/(2/i). The boundary 
conditions may be expanded in powers of i] and then e : 

4>0y\l + e[cOSa(x - t)ll> yy\l + ^ly\l] + g 2 ^ ^ 1 COS 2 Ol{x - t) + 4> 2 y\l + 

cosa(x - H = -Kn{^ 0ra |i + e [cos a (a; - t)i/j yyy\ l + + 

e 2 [^y^ cos 2 a (a; - 1) + cosa{x - + $2w|i] + •••}, (6) 

V'oxli + e[cos a(x - t)ip 0xy \i + ipi x \i] + e 2 ^^- cos 2 a(x - t) + 

cos a(x — t)tpixy\i + ip2x\i] + • • • = — ea sin a(x — i). (7) 
Equations above, together with the condition of symmetry and a uniform (dp/dx) , yield : 

^ = /To[(l + 2Kn)y-£], ^° = T ( ~^ )o ' (§) 

= \{^yy a{x - t] + 4> m (y)e-^*-% (9) 

where the asterisk denotes the complex conjugate. A substitution of ipi into Eqn. (3) yields 

d 2 d 2 
{— -a 2 + iaRell - K (l - y 2 + 2Kn)]}(— - a 2 )<p - 2iaK Re = 
dy 2 dy 2 

or if originally the (matter) gas is quiescent : K Q = (this corresponds to a free (vacuum) 
pumping case) 

(^-« 2 )(^7-« 2 )0 = O, a 2 = a 2 -iaRe. (10) 
dy 2 dy 2 

The boundary conditions are 

<f) y (±l) ± ra (±l)Kn = 2K (1 ± Kn) = 0, 0(±1) = ±1. (11) 

Similarly, with 

^2 = \{D{y) + E(y)e i2a ^ + E*(y) e - i2a ^} : (12) 



we have 

D yyyy = ~ 4>*<l>vy)v, (13) 

r/ 2 r/ 2 
[— - (4a 2 - 2Me)](— - Aa 2 )E - i2aReK (l - y 2 + 2Kn) 
ay 2 ay 2 

d 2 icxRc 
(jl ~ Aa2 ) E + ^aK ReE + -^-{^ yy - Hyyy) = 0; (14) 

and the boundary conditions 
D y (±l) + \[<i>yy(±l) + ^(±1)] - 2K = =^{±[^(±1) + <^(±1)] + £> w (±l)}, (15) 

E v (±l) + l<P yy (±l) - ^ = =FKn[^(±l) + ^(±1)], (16) 

E(±l) + i^(±l) = (17) 
where K is zero in Eqns. (13-16). After lengthy algebraic manipulations, we obtain 

= Co e ay + Cl e- ay + c 2 e" y + c 3 e-" y , 

where c = (A + A^/Det, c x = -(B + B )/Det, c 2 = (C + C )/Det, c 3 = — (T + T )/Det; 

Det = Ae a - Be~ a + Ce a - Te _ ", 



A = 


e a a 2 (r 2 e 


2d _ 


- s 2 e m ) - 2aae~ a w + aae a z{e~ 2& r + e m s), 


A) = 


e- a a 2 (r 2 e 


-2d 


- s 2 e 2a ) + 2aae a z - aae- a w(e m s + e" 2a r), 


B = 


e- a a 2 (r 2 e 


-2d 


- s 2 e 2a ) + 2aae a z - aae- a w(e~ 2a r + e 2a s), 


B = 


-- e a a 2 (r 2 e 


-2d 


- s 2 e m ) - 2aae~ a w + aae a 2(e" 2s r + e m s), 



C = e- a aa(wse*- a - rze a ~^) - ae 2a+a z(az - as) + ae- a w{ae«~ a w - ae a ~"r), 
C = e a aa(wse a - a - rze a ~ a ) - az(zae 2a ~ a - ae a s) + aw(ae-^ +2a) w - ae -(2a+a V), 
T = e - a aa(zse a+a - rae" (a+s) ) - aa(e 2a ~ & zr - e a ws) + a 2 e- a (-e 2a z 2 + e~ 2a w 2 ), 
T = e a aa(zse a+a - ™e" (a+s) ) - aa(e~°zr - e & ' 2a ws) + a 2 e a (-e 2a z 2 + e~ 2a w 2 ), 
with r = (1 — «Kn), s = (1 + aKn), w — (1 — «Kn), z = (1 + aKn). 

To obtain a simple solution which relates to the mean transport so long as only terms of 
0(e 2 ) are concerned, we see that if every term in the x- momentum equation is averaged over 
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an interval of time equal to the period of oscillation, we obtain for our solution as given by 
above equations the time-averaged (unit) body forcing 

I = ^ - 2 'tf + f«» ~ + °< e3 » = ^ + < 18 > 

where n is the integration constant for the integration of equation (12) and could be fixed 
indirectly in the coming equation (22). Now, from Eqn. (14), we have 

D y (±l) ± KnD yy (±l) = ~[4> w (±l) + <P* yy (±l)] T Kn{±[<f> yyy (±l) + <^(±1)]}, (19) 

where D y (y) = Hoy 2 + a±y + a 2 + C(y), and together from equation (12), we obtain 




with gi — a + a*, g 2 = a + a, g% = a — a*, g 4 = a — a, g 5 = a — a*, g 6 = a + a*. In 
realistic applications we must determine n from considerations of conditions at the ends of 
the matter-region. a\ equals to zero because of the symmetry of boundary conditions. 
Once n is specified, our solution for the mean speed (u averaged over time) of matter-flow 
is 

U = e 2 ^ = ^{C(y) - C(l) + i?o - KnC,(l) + U [y 2 - (1 + 2Kn)]} (21) 

where R = — {[0 W (1) + 0^(1)] _ Kn[0 TO2/ (l) +<f)* yy (l)]} /2, which has a numerical value 
about 3 for a wide range of a and Re (playing the role of viscous dissipations) when Kn= 0. 
To illustrate our results clearly, we adopt U(Y) = u(y) for the time-averaged results with 
y = Y in the following. 

III. RESULTS AND DISCUSSION 

We check our approach firstly by examining i?o with that of no-slip (Kn= 0) approach. 
This can be done easily once we consider terms of D y (y) and C(y) because to evaluate R 
we shall at most take into account the higher derivatives of 4>(y), like (j> yy (y), <Pyyyi.ll) instead 
of <j) y (y) and escape from the prescribing of a 2 . 
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Our numerical calculations confirm that the mean streamwise velocity distribution (averaged 
over time) due to the induced motion by the wavy elastic vacuum-matter interface in the 
case of free (vacuum) pumping is dominated by R (or Kn) and the parabolic distribution 
— n (l — y 2 ). Rq which defines the boundary value of D y has its origin in the y-gradient of 
the first-order streamwise velocity distribution, as can be seen in Eqn. (14). 
In addition to the terms mentioned above, there is a perturbation term which varies across 
the channel : C(y) — C(l). Let us define it to be 

F{y) = ^ [C{y) - C{l)] (22) 

To compare with no-slip (Kn= 0) results, we plot three cases, a = 0.1,0.4, and 0.8 for the 
same Reynolds number Re = 1 of our results : Kn= 0.1 with those Kn= into Fig. 1. We 
remind the readers that the Reynolds number here is based on the wave speed. This figure 
confirms our approaches since we can recover no-slip results by checking curves of Kn= 
and finding them being almost completely matched in-between. The physical trend herein 
is also the same as those reported in Refs. [12-13] for the slip-flow effects. The slip produces 
decoupling with the inertia of the wavy interface. 

Now, let us define a critical reflux condition as one for which the mean velocity U (Y) equals 
to zero at the center-line Y = (cf. Fig. 1). With equations (12,20-21), we have 



= fi A = [^F(0)/2QQ + Kn C -(l)-^,l 
0cr W -(l + 2Kn) V ; 

which means the critical reflux condition is reached when ITo has above value. Pumping 
against a positive (unit) body forcing greater than the critical value would result in a back- 
ward transport (reflux) in the central region of the stream. This critical value depends on 
a, Re, and Kn. There will be no reflux if the (unit) body forcing or pressure gradient is 
smaller than this n . Thus, for some n values less than n 0cr , the matter (flow) will keep 
moving or evolving forward. On the contrary, parts of the matter (flow) will move or evolve 
backward if Ho > Ho cr . This result could be similar to that in [16] using different approach 
or qualitatively related to that of [5] : even for very slow growth of A, the gravitationally 
bound systems become unbound while the nongravitationally bound systems remain bound 
for certain parameters defined in [5] (e.g., rf). 

We present some of the values of n (a,-Re;Kn = 0,0.1) corresponding ot freezed or zero- 
volume-flow-rate states (j\U(Y)dY = 0) in Table 1 where the wave number (a) has the 



range between 0.20 and 0.80; the Reynolds number (Re)— 0.1, 1, 10, 100. We observe that as 
Kn increases from zero to 0.1, the critical n or time-averaged (unit) body forcing decreases 
significantly. For the same Kn, once Re is larger than 10, critical reflux values Il drop 
rapidly and the wave-modulation effect (due to a) appears. The latter observation might 
be interpreted as the strong coupling between the vacuum-matter interface and the inertia 
of the streaming matter-flow. The illustration of the velocity fields for those zero-flux (zero- 
volume-flow-rate) or freezed states are shown in Figure 2. There are three wave numbers 
: a = 0.2,0.5,0.8. The Reynolds number is 10. Both no-slip and slip (Kn=0.1) cases are 
presented. The arrows for slip cases are schematic and represent the direction of positive 
and negative velocity fields. 

Some remarks could be made about these states : the matter or universe being freezed 
in the time-averaged sense for specific dissipations (in terms of Reynolds number which is 
the ratio of wave-inertia and viscous effects) and wave numbers (due to the wavy vacuum- 
interface or vacuum fluctuations) for either no-slip and slip cases. This particular result 
might also be related to a changing cosmological term (growing or decaying slowly) or the 
critical density mentioned in Refs. [1-2]. If we treat the (unit) body forcing as the pres- 
sure gradient, then for the same transport direction (say, positive x-direction), the negative 
pressure (either downsdtream or upstream) will, at least, occur once the time-averaged flow 
(the maximum speed of the matter (gas) appears at the center-line) is moving forward! 
For example, for n = Re(^) 2 = —10 < (Re — l,a — 0.5), the velocity field (pro- 
file) is shown in Fig. 3. One possible p-pair for uniform (negative) gradient (mean value 

theorem). Pdownstream Pupstream ^ 0, with Xdownstream Xupstream ^ . Pdownstream *^ 0, 
Pdownstream *^ Pupstream *^ 0. 

Meanwhile, the time-averaged transport induced by the wavy interface is proportional to 
the square of the amplitude ratio (although the small amplitude waves being presumed), 
as can be seen in Eqn. (12) or (20), which is qualitatively the same as that presented in 
[9] for analogous interfacial problems. In brief summary, the entrained transport (pattern, 
either postive or negative and there is possibility : freezing) due to the wavy vacuum-matter 
interface is mainly tuned by the (unit) body forcing or n for fixed Re. Meanwhile, n 0cr 
depends strongly on the Knudsen number (Kn, a rarefaction measure) instead of Re or a. 
We hope that in the future we can investigate other issues like the role of phase transition 
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and that of cyclic universes [16-18] using the present or more advanced approach. 
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U(Y) ! 




Fig. 1 Demonstration of Kn and a effects on the time-averaged velocity (fields) profiles. 
IIo = no cr . The Reynolds number (the ratio of the wave- inertia and viscous dissipation 
effects) is 1. U(Y) = at Y = 0. Kn is the Knudsen number (a rarefaction measure). 
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Fig. 2 Demonstration of the zero- flux states : the mean velocity field UiY) for 
wave numbers a = 0.2, 0.5, 0.8. The Reynolds number is 10. Kn is the rarefaction measure 
(the mean free path of the particles divided by the characteristic length) . 
The arrows are schematic and illustrate the directions of positive and negative U(Y). 
The integration of U(Y) w.r.t. Y for these velocity fields gives zero volume flow rate. 
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Re=1.0, n Q =-10; o : Kn=0, + Kn=0.1 o 
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Fig. 3 Demonstration of the negative-p states : the mean velocity field U (Y) 
for the wave number a = 0.5 and the Reynolds number Re = 1. 
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TABLE I: Zero-flux or freezed states values (IIq) for a flat vacuum-matter interface. 







Re 








Kn 


a 


0.1 


1 


10 


100 





0.2 


4.5269 


4.5269 


4.5231 


4.3275 




0.5 


4.6586 


4.6584 


4.6359 


4.2682 




0.8 


4.9238 


4.9234 


4.8708 


4.4488 


0.1 


0.2 


2.4003 


2.4000 


2.3774 


1.2217 




0.5 


2.4149 


2.4132 


2.2731 


-0.9054 




0.8 


2.4422 


2.4379 


2.0718 


-3.4151 
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